Phase diagrams of period-4 spin chains consisting of three kinds of spins 
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We study a period-4 antiferromagnetic mixed quantum spin chain consisting of three kinds of spins. 
When the ground state is singlet, the spin magnitudes in a unit cell are arrayed as (s — t, s, s + 1, s) 
with integer or half-odd integer s and t (0 < t < s). The spin Hamiltonian is mapped onto a 
nonlinear a model (NLSM) in a previously developed method. The resultant NLSM includes only two 
independent parameters originating from four exchange constants for fixed s and t. The topological 
angle in the NLSM determines the gapless phase boundaries between disordered phases in the 
parameter space. The phase diagrams for various s and t shows rich structures. We systematically 
explain the phases in the singlet-cluster-solid picture. 



PACS numbers: 75.10.Jm, 75.30. Kz 
I. INTRODUCTION 

The lowest spin excitation for a homogeneous antifer- 
romagnetic spin chain has a gap if the spin magnitude is 
an integer, and has no gap if it is a half odd integer. Hal- 
dane conjectured this proposition by mapping the spin 
chain onto the nonlinear a model (NLSM) p| . From the 
topological term of the NLSM, we know whether the spin 
excitation is gapless or not. Mapping to the NLSM is now 
recognized as one of powerful methods to examine quan- 
tum spin systems. Various aspects of the NLSM method 
are found in Refs. §,||]. 

Analyses using the NLSM have been extended to in- 
homogeneous spin chains. The first step has been done 
by Affleck ffl. He reformulated the NLSM method in 
an operator formalism in which spin operators are trans- 
formed in pairs. Due to the pair transformation, his for- 
malism is applicable to a spin chain with bond alterna- 
tion. Then the exchange constant is inhomogeneous with 
period 2 but the spin magnitude is still homogeneous. 
Results by the NLSM method qualitatively agree with 
numerical and experimental results fj],|],f§ ■ 

The NLSM method for a general inhomogeneous spin 
chain with arbitrary period is given, if the ground state 
is singlet, in a previous paper |10|. The inhomogeneity is 
not only for the exchange constant but also for the spin 
magnitude. The derivation of the NLSM is based on 
dividing the spin chain into blocks and simultaneously 
transforming the spin variables belonging to a block in a 
path integral formalism. From the topological term, the 
gapless condition for the spin excitation is obtained as 
an equation, which we have called the gapless equation. 
The gapless equation determines the phase boundaries 
in the ground-state phase diagram. This method is quite 
general, and is applicable to various spin chains. 

A period-4 chain is the simplest nontrivial case in 
which more than one kind of spins can be mixed un- 
der condition that the ground state is singlet. We have 
applied the NLSM method to a period-4 chain consisting 
of two kinds of spins and presented the phase diagram in 
the parameter space of exchange constants . When 



the magnitudes of two kinds of spins are s a and Sf,, only 
the array of (s a , s a , s&, Sb) in a unit cell consists with 
a singlet ground state. We have constructed an NLSM 
and then the gapless equation for this case following the 
general procedure flfj|| ; Fukui and Kawakami also derived 
the NLSM in a different standpoint |l3| . We constructed 
phase diagrams from the NLSM and found that they gen- 
erally consist of many disordered phases. To understand 
various phases we have proposed the singlet-cluster-solid 
(SCS) picture, which is an extension of the valence-bond- 
solid (VBS) picture |l4"| . The SCS picture systematically 
explains all phases for any values of s a and s&. A simple 
version of the SCS picture is seen in the case of s a =Sb=h 
p"5| . In the case of s a =iand Sb=l, numerical calcula- 
tion has been performed jl6| ] . The phase diagram by the 
NLSM method qualitatively agrees with the numerical 
phase diagram. 

In this paper, we study a period-4 spin chain consist- 
ing of three kinds of spins. Systematic treatment of this 
problem seems to be difficult due to many possibilities 
in the choice of the spin kinds. However, the condition 
that the ground state is singlet fairly restricts allowed 
combinations. In fact, the spin magnitudes in a unit cell 
must be generally arrayed as (s — t, s, s + t, s) with in- 
teger or half-odd integer s and t (0 < t < s), as will 
be seen. The spin chain has a modulation in the spin 
magnitude around average s; hence it is interesting to 
compare results for different t and the same s. The ex- 
change constants are also periodic with period 4 and then 
there are three parameters except for an energy unit. 
The spin Hamiltonian is mapped onto an NLSM in the 
general method jl0|]. The resultant NLSM shows that 
the number of relevant parameters is not three but two. 
The phase diagrams for various s and t determined by 
the gapless equation show rich phase structures (ijj . We 
systematically explain the phases in the SCS picture. 

This paper is organized as follows. In Sec. II, we in- 
troduce the Hamiltonian for the period-4 spin chain con- 
sisting of three kinds of spins with a singlet ground state, 
and parameterize the exchange constants. In Sec. Ill, the 
spin Hamiltonian is transformed to an NLSM. From the 
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topological term of the NLSM, the gapless equation to 
determine phase boundaries is derived. In Sec. IV, phase 
diagrams are drawn by the gapless equation for various 
spin magnitudes. Features of the phases are mentioned. 
In Sec. V, the ground states of the phases are explained 
by the SCS picture. Section VI is devoted to summary 
and discussion. 
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II. HAMILTONIAN 

The spin Hamiltonian with period 4 is generally writ- 
ten as 

N/4 

S47+I ■ S4j + 2 + J2 S4J+2 ' S4J+3 



J4 S 4j+ 4 ■ S 4j+5 ), (1) 



3 °4j+3 ' >->4j + 4 ~r J 4 &4j+4 1 a 4j+5 

where is the spin at site j with magnitude Sj. The 
number of lattice sites is N, the lattice spacing is a and 
the system size is L = aN. We only consider the antifer- 
romagnetic exchange interaction (J,; > 0). 

Since the system is periodic with period 4 in the spin 
magnitude as well as in the exchange constant, we gen- 
erally have 4 different spin magnitudes (s 4 , S2, S3, S4) in 
a unit cell. A unit cell is illustrated in Fig. The values 
of the spin magnitudes are restricted in order that the 
system has a singlet ground state. Following the Lieb- 
Mattis theorem p8| , a singlet ground state is realized if 
the system satisfies the restriction 



si - s 2 + s 3 - S4 = 0; 



(2) 



otherwise the system has a ferrimagnetic ground state. 
We have treated the case that three kinds of spins are 
mixed. To consist with the restriction (0), the spin mag- 
nitudes must be of the following form: 



(si, s 2 , s 3 , s 4 ) = (s - t,s,s + t, s), 



(3) 



where s and t are positive integers or half-odd integers 
satisfying s > t > 0. The spin configuration (^) is re- 
garded as a modulation against the uniform configuration 
(s,s,s,s). 

We parameterize the exchange constants as 



J 1 = 



J 



1 + 5' 
J' 



J 2 = 



J' 



J3 — ^— — j J4 — 
1 + 7 



1-7 
J 

1~ 7' 



(4) 



where S (7) is the distortion parameter describing the 
asymmetry between exchange constants J4 and J± (J2 
and J3) on the both sides of a spin with s± (S3). 



FIG. 1. Spin magnitudes and exchange constants of the 
Hamiltonian (1) with Eqs. (3) and (4) in a unit cell. 



III. MAPPING TO THE NONLINEAR a MODEL 

The expectation value of a spin operator in a coherent 
state is expressed as 



< S, >= (-ly'ajiij 



(5) 



with a unit vector rij . The partition function Z at tem- 
perature 1//3 is then represented as 



(6) 



Z= jD[n 3 ] J]*(n?-l)e- s , 

j 

N 

1 [ p N 

+ 2/ dr IJ J J s i s J+i( n J - n J+i) 2 - ( 7 ) 



In the action (0), the first term comes from the Berry 
phase and uj[nj] is the solid angle which the unit vector 
iij forms in the period (3. 

We have derived the NLSM action starting from the 
action (^) for the general periodic case pp| . The deriva- 
tion is based on dividing the spin chain into blocks and 
transforming the spin variables in a block into new ones. 
In the present case, by choosing unit cells as blocks, the 
transformation for the pth block is written as 

3 1 
n 4p +i = jm(p) + -m(> - 1) + oLi(p), 



m(p) 



•oL 2 (p), 



3 1 
n 4p +3 = |m(p) + -m(p + 1) + aL 3 (p), 

n 4p +4 = 2 m (P) + 2 m ( p + X ) + aL 4b)> (8) 

where {m(p)} are gradually changing unit vectors and 
{L 9 (p)} are small fluctuations. This transformation does 
not change the number of the original degrees of freedom 
p0| . Integrating out the fluctuations {L q (p)} and taking 
the continuum limit, we obtain the effective action p9|: 
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rfi rL f j(0) 

S C B = J dr J dxl-i-j^ra ■ (d T m x d x m) 

^(8 r m) a + |j(°)(fltm) a }. 



J(0) 



1 



2aJW V J( 2 ) 
where 
1 1 

1 _ s - t 
JW ~ IT 

1 1 
jW _ 4s(s + t) 



1 



J(s-i) J'(s + t) 
1 1 



J(s-t) J'(s + t) J'(s-t) 



s 2 + f 



J J' J' 

We have used the effective distortion parameter 

r J ' 

d = j + S—. 

It is remarkable that two distortion parameters 7 and S 
appear only in a single parameter d. Hence the system 
is characterized by only a pair of parameters, (J'/J, d). 
The value of d is restricted as 



(9) 

(10) 

(11) 
(12) 

(13) 



\d\ <1 



(14) 



because \5\ < 1 and I7I < 1 from Eq. (||). 

The action (O) is of the standard form of the NLSM: 



■P rL 

dr I dx 
Jo 



-i — m ■ (d T m x d x m) 
Air 



(15) 



with the topological angle 8, the coupling constant g and 
the spin wave velocity v. Comparing Eq. (^) to Eq. (15), 
9 is given as 



1 = 4tt 



J« 



(16) 



IV. PHASE DIAGRAMS 



The NLSM has a gapless excitation when 0/2tt is a 
half-odd integer. From Eq. (16), this condition is rewrit- 
ten as the gapless equation: 



1 _2/o-l 1 
JW ~ _ 



J(°) 



(17) 



where lo is an integer. For each lo, this equation deter- 
mines a boundary between disordered phases if it has a 
solution. 



Substituting Eqs. (|o|) and ([TT]) into Eq. (p^|), we have 
1\ / 1 1 J r 



I 



t 



s + t s — t J 



(18) 
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FIG. 2. The phase diagrams of s = 1 for (a) t = 0, and (b) 
|. Unshaded regions are physical. 

with i = lo + 2t + 1. A phase boundary is a straight line 
in the parameter space of (J'/J, d). Allowed values of 
integer I satisfying Eq. ([14]) are 



I = 1,2, 



,2(s + t). 



(19) 



Thus we have found that there are 2(s + 1) + 1 phases 
separated by 2(s + t) gapless phase boundaries. All of 
them go through the point 



s - t 
' s + f 



,0 



(20) 



In the special case of d = 0, Eq. ( Jlq ) reduces to Z = 
s + £ + h. Then, the system has a spin gap when S3 
(= s + t) is an integer, while it has no spin gap when 
S3 is a half-odd integer. This is an extended version of 
Haldane's proposition for t = and J' = J to cases with 
period-4 modulation with respect to the spin magnitude 
and the exchange interaction. The condition d = (i. e. 
J7 + J'S = 0) is satisfied even for nonzero 7 and 5. This 
means that effects of distortions can be canceled. 

In Fig. ^, we present phase diagrams for 8 = 1. Figure 
||(a) is for s = 1 and t = 0; spins are arrayed as 1-1-1-1 in a 
unit cell. Then the spin magnitude is homogeneous, while 
the exchange interaction is modulated with period 4. The 
shaded regions are not physical because of Eq. ([l4]) and 
J'/J > 0. There are three phases labeled by A + , H 7 
and A- . Phase H is the Haldane phase since it includes 
the point (1,0) which represents the uniform s = 1 spin 
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FIG. 3. The phase diagrams of s = | for (a) t = 0, (b) |, 
and (c) 1. Unshaded regions are physical. 



chain. Two phases A + and A_ appear, when the effective 
distortion is strong. Figure ||(b) is for s = 1 and t = 
^; spins are arrayed as in a unit cell. There are 

four phases labeled by B + , A + , A_ , and £?_ . In contrast 
to the case of (a), the system is on a boundary when 
there is no effective distortion (d=0). It is noticed that 
the phases B + and J3_ appear, when J'/J < 1. All the 
phases in Figs. ||(a) and (b) are interpreted by means of 
the SCS picture in the next section. 

In Fig. ||, we present phase diagrams for s=|. Figure 
||(a) is for s=| and t=0; spins are arrayed as §-§-§-§ in 
a unit cell. Only the exchange interaction is modulated 
with period 4. There are 4 phases labeled by B + , A + , 
A_, and 73_. These phases, including B + and £?_ for 
strong distortion, develop from small J'/J to large J'/J. 
Figure ||(b) is for s=| and t=\\ spins are arrayed as 
1-I-2-I in a unit cell. There are 5 phases labeled by B + , 
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FIG. 4. The phase diagrams of s = 2 for (a) i = 0, (b) |, 
(c) 1, and (d) |. Unshaded regions are physical. 
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A+, H, A. 



and B_. Phases B + and B_ appear only 
when J'/J < 1. Figure ||(c) is for s=| and t=l; spins 
are arrayed as |~f~§~| in a unit cell. There are 6 phases 
labeled by C+, B +1 A + , A-, B-, and C_. Phases C+, 
£? + , and C_ appear only when J'/J < 1; C+ and 
C_ are for strong distortion and for very small J' / J. In 
the case of (a) , the system is on a gapless boundary when 
there is no effective distortion (d=0). All the phases in 
Figs. ||(a), (b) and (c) are interpreted by means of the 
SCS picture in the next section. 

In Fig. |[ we present phase diagrams for s=2. Figures 
4(a), (b), (c) and (d) are for t—0, t=h, t=l, and i=§, 
respectively; spins are arrayed as 2-2-2-2, |-2-|-2, 1-2- 
3-2, and |-2-|-2 in a unit cell. Phase H for small d 
in Figs. 4(a) and (c) is the s—2 Haldane phase, which 
includes the point (1,0). In the cases of (b) and (d), the 
system is on a gapless boundary when there is no effective 
distortion (g?=0). Some phases develop from small J'/J 
to large J'/ J, and the others are restricted in the region 
J'/J < 1. 



V. SINGLET-CLUSTER-SOLID PICTURE 

Disordered phases in the phase diagrams are explained 
in the SCS picture, which includes the VBS picture as a 
special case. In the SCS picture, a spin with more than 
1 magnitude is decomposed into h spins. The original 
spin state is retrieved by symmetrizing the states of the 
decomposed h spins at each site. A disordered state is 
represented by a regular array of local singlet clusters of 
even numbers of \ spins, while a state in the VBS picture 
is a direct product of only singlet dimers. 

The SCS pictures for the phase diagrams of s=l in 
Fig. U are shown in Figs, 
the phase diagram in Fig. I 

(s=l, t=0). For each picture, a small circle represents a \ 
spin; an original s—1 spin is decomposed into two i spins. 
A loop represents a singlet dimer of two \ spins in it. A 
dashed line represents a spatially extended singlet state. 
Picture (a) ((e)) is for a state in the dimer phase A + 
(yl_), and has advantage for the exchange energy on J2 
and/or J4 (Ji and/or J 3 ) interactions because of positive 
(negative) d. Picture (c) is for a translationally invariant 
state in the Haldane phase H . These pictures are nothing 
but the VBS picture; a ground state is represented only 
by singlet dimers. Pictures (b) and (d) are for states 
between A + and H , and between H and A- , respectively; 
each picture includes an extended singlet state (dashed 
line) contributing to a gapless excitation. 

Figure |^ explains the phase diagram in Fig. ^|(b) for 
the |-1-|-1 spin chain (s=l, t—^). Pictures (c) and (e) 
are for dimer states A + and A- , respectively. The state 
(c) ((e)) has advantage for the exchange energy on J 2 
and/or J4 (Ji and/or J3) interactions because of posi- 
tive (negative) distortion d. Pictures (a) and (g) are for 
singlet cluster states in the phases B + and £>_, which 



Ji J2 J3 J4 



and ^. Figure |5| explains 
2(a) for the 1-1-1-1 spin chain 



(a) 



(b) 



3) 


Co 


O) 


Co 


0) 


Co 


0) 


en 


Zo) 


(0 


O) 


(0 


0) 


(0 


0) 


CqI 


-0- 


--0-- 


-O- 


-0- 


--0- 


-0-- 


-0- 


-0- 


Zo) 


Co 


O) 


(0 


0) 


(0 


0) 


(si 


Co 


0) 


Co 


0) 


Co 


0) 


Co 


0) 


Zo) 


(n 
\P 


O) 






(P 


V) 




Co 


O) 


(0 


0) 


(0 


0) 


Co 


0) 


-0- 


--0-- 


-0- 


-0- 


- -CD- 


-0-- 


-0- 


-0- 


1 PI 


O) 


Co 


0) 


CO 


0) 


Co 


0) 


(0 


O) 


(0 


0) 


(0 


0) 


(0 


0) 



FIG. 5. The SCS pictures for the phases of s=l and t~0. 
They represent states for phases (a) A+, (c) H, and (e) A- 
in Fig. 2(a). H is the Haldane phase. They are the same as 
the VBS pictures in the present case. Gapless states on the 
phase boundaries are presented in (b), and (d). 



are in the region of J'/J < 1. In phase B + (-B-), a sin- 
glet cluster or dimer including a J3 (J2) interaction is 
the most unfavorable for the exchange energy because of 
large \d\ and small J'/J. Hence clusters of six 1 spins 
are formed to avoid singlet dimers at J3 (J2) interac- 
tions. Pictures (b), (d), and (f) are for gapless states on 
the phase boundaries, where extended states appear. 

The SCS pictures for the phase diagrams of s=| in 
Fig. U are shown in Figs. [?], ||, and ^. Those for phase 
boundaries are not drawn to reduce the figure sizes. Fig- 
ure [7] explains the phase diagram in Fig. |3](a) for the 
s P^ n chain (s=|, t—0). The SCS pictures repre- 
sent the ground states in phases (a) B + , (b) A + , (c) ^4_, 
and (d) B_, and are the same as the VBS pictures. Pic- 
ture (a) ((d)) is for the dimer phase B + (B_), where the 
energy reduction on J2 and/or J4 (Ji and/or J3) interac- 
tions is the most favorable because of large \d\. Picture 
(b) ((c)) is for the phase A + (A_), where the energy re- 
duction on J2 and/or J4 (Ji and/or J3) interactions is 
favorable to some extent because of smaller but finite \d\. 

Figure explains the phase diagram in Fig. ||(b) for 
the 1-|-2^| spin chain (s=|, t—^). The SCS pictures 
represent states in phases (a) B + , (b) A + , (c) H, (d) A_, 
and (e) Picture (c) represents a Haldane-like state. 
With transitions (c) to (b) and (b) to (a), the number of 
dimers including J3 interactions decreases. This explains 
that d in B + is larger than d in A + , and d in A + is larger 
than d in H as seen in Fig. ||(b). In particular, the SCS 
pictures in B + and £?_ include singlet clusters of six i 
spins. The large clusters are formed for J'/ J < 1 to avoid 
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FIG. 6. The SCS pictures for the phases of s=l and t=g. 
They represent states for phases (a) j3_|_, (c) t4 + , (e) _4_, and 
(g) -B- in Fig. 2(b). Gapless states on the phase boundaries 
are presented in (b), (d), and (f). 



(b) 



Co o) (o o) Co o) Co o) 
o) Col 
o) CH 



~^6) (o o) (o o) 
^) (o o) (o o) 



(o o) Co ci) (o o) (o o) 



Co o) (o o) (o o) (o o) 



(o o) (o o) (o o) 



(d) (o o) (o o) (o o) (o o) 



(o o) (o o) (o o) (o o) 



(o o) (o o) Co o) Co o) 

FIG. 7. The SCS pictures for the phases of s=| and t=0. 
They represent states for phases (a) B+, (b) A+, (c) j4_, and 
(d) B_ in Fig. 3(a). They are the same as the VBS pictures 
in the present case. 



singlet dimers including J3 (J2) interactions. 

Figure explains the phase diagram in Fig. ||(a) for 
the |-|-|-| spin chain (s=|, t=l). The SCS pictures 
represent states in phases (a) C + , (b) B + , (c) A + , (d) 
_4_, (c) j5_, and (f) C_. Large singlet clusters of more 
than two | spins are formed in C + , S + , and C_ for 
J'/ J < 1 again. 

The SCS pictures of the ground states for arbitrary 
s and t, including the cases of s=2 in Fig. [|, are sim- 
ilar to the above examples. Phases for small values of 
\d\ are represented by VBS pictures, each of which con- 
sists of singlet dimers only. The number of possible VBS 
pictures is 2(s — t), and they are between the bound- 
aries of l=2t and of l=2s + 1 in Eqs. (|l|) and (|l|). For 
larger d outside the VBS phases in the (J' /J, d) space, 
there appear phases represented by SCS pictures which 
includes singlet clusters larger than dimers. Generally, 
as one moves from a phase to another accompanied by 
increasing (decreasing) d for d > (d < 0), the num- 
ber of singlet dimers on J3 ( J2) interactions in the SCS 
picture decreases by one per unit cell. This explains the 
total number 2(s + t) + 1 of the phases. The phases 
explained by SCS pictures including singlet clusters con- 
sisting of more than two ^ spins are restricted to small 
J' I J regions. In fact, singlet dimers including J 3 (J 2 ) 
interactions are energetically unfavorable for small J'/ J, 
and the number of them can be reduced for t ^ if large 
clusters not including J3 (J2) interactions are formed. 
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FIG. 8. The SCS pictures for the phases of s=| and t=\. 
They represent states for phases (a) B + , (b) A + , (c) H, (d) 
A-, and (e) B- in Fig. 3(b). 
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FIG. 9. The SCS pictures for the phases of s=§ and t=l. 
They represent states for phases (a) C+, (b) B + , (c) A + , (d) 
A-, (e) B+, and (f) C_ in Fig. 3(c). 
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VI. SUMMARY AND DISCUSSION 



We studied a period-4 antiferromagnetic mixed quan- 
tum spin chain consisting of three kinds of spins in the 
case that the ground state is singlet. The spin magni- 
tudes in a unit cell must be arrayed as (s — t, s, s + t, s) 
with integer or half-odd integer s and t (0 < t < s). 
The exchange constants in a unit cell are (Ji, J2, J3, J4); 
i. e. there are three parameters except for an energy 
unit. The spin Hamiltonian is transformed to an NLSM 
by a previously developed method [0. In the NLSM, 
we find that the number of relevant parameters is not 
three but two. One of them is a ratio J' / J of the ex- 
change energies without distortion, and the other is the 
effective distortion parameter d. By the gapless equation 
from the NLSM, we determined 2s + 2t phase boundaries 
between gapful disordered phases in the (J' /J, d) space 
for each pair of s and t. We explained systematically the 
disordered phases by means of the SCS pictures. In the 
case of d = 0, the ground state is in a gapful Haldane- 
like phase for s + t being an integer, and is on a gapless 
phase boundary for s + t being a half-odd integer. When 
d increases (decreases) from d — 0, singlet clusters in- 
cluding J 3 (J 2 ) interactions successively decreases with 
phase transitions. If t ^ 0, singlet clusters larger than 
dimers appear for large |d| and small J' /J. 

We discuss whether a change from an SCS state to 
another is a phase transition or not. For simplicity, we 
consider a part of a full SCS state as shown in Fig. [II]. 
State (a) is a VBS state which is a direct product of sin- 
glet dimers. This state can gradually change to state (b) 
without a phase transition, because state (b) is formed 
by local modifications where a pair of dimers puts to- 
gether into a singlet cluster of four spins. The change 
from (a) to (c) is similar. We have used picture (a) as a 
representative of (a), (b) and (c) in this paper. On the 
other hand, there is no way to locally modify dimers in 
(a) to form dimers in (d). That is, the change from (a) 
to (d) must be realized only by a global recombination of 
dimers or a phase transition. The wave function for (a) is 
symmetric but that for (d) is antisymmetric with respect 
to the spatial reflection about the vertical dotted line. 
An extended state in (e) appears under the transition 
where both the dimer states are collapsed. 

Finally it is expected that materials realizing period- 
4 quantum spin chains will be synthesized and experi- 
mentally studied. The present paper (and Ref. [^2)) will 
hopefully work as a guide to investigate such materials. 
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